{smcl}
{com}{sf}{ul off}{txt}{.-}
      name:  {res}<unnamed>
       {txt}log:  {res}/N/project/suicide_study/pnas_replication/results/log/mi_1.smcl
  {txt}log type:  {res}smcl
 {txt}opened on:  {res}20 Aug 2020, 04:47:21
{txt}
{com}. 
. if ("`model'" == "logit"){c -(}
.         use "${c -(}home_dir{c )-}/data/processed/suicide_reg_v1_raw.dta", clear
. {c )-}
{txt}
{com}. 
. if ("`model'" == "mi") {c -(}
.         use "${c -(}home_dir{c )-}/data/processed/suicide_reg_v1_imputed_M10.dta", clear  
. {c )-}
{txt}
{com}. 
. * for now, we use the following simple survey weights 
. if ("`model'" == "mi") {c -(}
.         mi svyset `geo_type' [pw=ObsWgt0] 
{res}
      {txt}pweight:{col 16}{res}ObsWgt0
          {txt}VCE:{col 16}{res}linearized
  {txt}Single unit:{col 16}{res}missing
     {txt}Strata 1:{col 16}<one>
         SU 1:{col 16}{res}county
        {txt}FPC 1:{col 16}<zero>
{p2colreset}{...}
{res}{com}. {c )-}
{txt}
{com}. else {c -(}
.         svyset `geo_type' [pw=ObsWgt0]  
. {c )-}
{txt}
{com}. 
. * margins for each category
. program margin_interact 
{txt}  1{com}.         args X Y k model
{txt}  2{com}.         sum `X', d 
{txt}  3{com}.         local gap = (`r(max)' - `r(min)') / `k' 
{txt}  4{com}.         if ("`model'" == "logit") {c -(}
{txt}  5{com}.                 margin `Y', at(`X' = (`r(min)' (`gap') `r(max)')) predict(pr)
{txt}  6{com}.         {c )-} 
{txt}  7{com}.         else if ("`model'" == "mi") {c -(}
{txt}  8{com}.                 mimrgns `Y', at(`X' = (`r(min)' (`gap') `r(max)')) predict(pr)
{txt}  9{com}.         {c )-}       
{txt} 10{com}. end 
{txt}
{com}. 
. program mchange_mi
{txt}  1{com}.         args X k model
{txt}  2{com}.         if ("`model'" == "logit") {c -(}
{txt}  3{com}. 
.                 if ("`k'" == "continuous") {c -(}
{txt}  4{com}.                         sum `X' if e(sample), d 
{txt}  5{com}.                         margin, at(`X' = (`r(min)' `r(max)')) post predict(pr)
{txt}  6{com}.                         mlincom  2 - 1, decimal(7) stat(all)    
{txt}  7{com}.                 {c )-} 
{txt}  8{com}.                 else if ("`k'" == "binary") {c -(}
{txt}  9{com}.                         sum `X' if e(sample), d 
{txt} 10{com}.                         margin, at(`X' = (`r(min)' `r(max)')) post predict(pr)
{txt} 11{com}.                         mlincom  2 - 1, decimal(7) stat(all)
{txt} 12{com}.                 {c )-} 
{txt} 13{com}.                 else if ("`k'" == "categorical") {c -(}
{txt} 14{com}.                         margin `X' if e(sample)==1 , at() pwcompare predict(pr) 
{txt} 15{com}.                 {c )-}
{txt} 16{com}.         {c )-}
{txt} 17{com}. 
.         else if ("`model'" == "mi") {c -(}
{txt} 18{com}. 
.                 if ("`k'" == "continuous") {c -(}
{txt} 19{com}.                         sum `X' , d 
{txt} 20{com}.                         mimrgns, at(`X' = (`r(min)' `r(max)')) post predict(pr)
{txt} 21{com}.                         mlincom  2 - 1, decimal(7) stat(all)
{txt} 22{com}.                 {c )-} 
{txt} 23{com}.                 else if ("`k'" == "binary") {c -(}
{txt} 24{com}.                         sum `X' , d 
{txt} 25{com}.                         mimrgns, at(`X' = (`r(min)' `r(max)')) post predict(pr)
{txt} 26{com}.                         mlincom  2 - 1, decimal(7) stat(all)
{txt} 27{com}.                 {c )-} 
{txt} 28{com}.                 else if ("`k'" == "categorical") {c -(}
{txt} 29{com}.                         mimrgns `X' , at()   pwcompare predict(pr)      
{txt} 30{com}.                 {c )-}
{txt} 31{com}.         {c )-}
{txt} 32{com}. end
{txt}
{com}. 
. * create some dummy codings
. tab Race5, gen(race5_nh)

      {txt}Race5 {c |}      Freq.     Percent        Cum.
{hline 12}{c +}{hline 35}
          1 {c |}{res}107,831,845       79.53       79.53
{txt}          2 {c |}{res} 14,631,144       10.79       90.32
{txt}          3 {c |}{res}  2,867,304        2.11       92.44
{txt}          4 {c |}{res}  2,753,311        2.03       94.47
{txt}          5 {c |}{res}  7,499,129        5.53      100.00
{txt}{hline 12}{c +}{hline 35}
      Total {c |}{res}135,582,733      100.00
{txt}
{com}.         rename race5_nh1 White_nh 
{res}{txt}
{com}.         rename race5_nh2 Black_nh 
{res}{txt}
{com}.         rename race5_nh3 AIAN_nh 
{res}{txt}
{com}.         rename race5_nh4 AsPI_nh 
{res}{txt}
{com}.         rename race5_nh5 Hispanic
{res}{txt}
{com}. 
. tab MarStat5, gen(ms)

   {txt}MarStat5 {c |}      Freq.     Percent        Cum.
{hline 12}{c +}{hline 35}
          1 {c |}{res} 75,689,936       55.83       55.83
{txt}          2 {c |}{res} 10,166,664        7.50       63.32
{txt}          3 {c |}{res} 13,925,509       10.27       73.60
{txt}          4 {c |}{res}  2,609,308        1.92       75.52
{txt}          5 {c |}{res} 33,190,855       24.48      100.00
{txt}{hline 12}{c +}{hline 35}
      Total {c |}{res}135,582,272      100.00
{txt}
{com}.         rename ms1 Marrd5
{res}{txt}
{com}.         rename ms2 Widow5
{res}{txt}
{com}.         rename ms3 Divor5
{res}{txt}
{com}.         rename ms4 Separ5
{res}{txt}
{com}.         rename ms5 NvMar5
{res}{txt}
{com}. 
. tab AgeGrp4, gen(ag) 

    {txt}AgeGrp4 {c |}      Freq.     Percent        Cum.
{hline 12}{c +}{hline 35}
          1 {c |}{res} 21,231,562       15.66       15.66
{txt}          2 {c |}{res} 39,043,312       28.80       44.46
{txt}          3 {c |}{res} 47,351,260       34.92       79.38
{txt}          4 {c |}{res} 27,956,599       20.62      100.00
{txt}{hline 12}{c +}{hline 35}
      Total {c |}{res}135,582,733      100.00
{txt}
{com}.         rename ag1 Age_15_24
{res}{txt}
{com}.         rename ag2 Age_25_44
{res}{txt}
{com}.         rename ag3 Age_45_64
{res}{txt}
{com}.         rename ag4 Age_65_Up
{res}{txt}
{com}. 
. destring St, replace 
{txt}St: all characters numeric; {res}replaced {txt}as {res}byte
{txt}
{com}. 
. * set-up equations
. local religion Rat_GC_ProE Rat_GC_ProM Rat_GC_ProB Rat_GC_Cath Rat_GC_Jew Rat_GC_Oth
{txt}
{com}. local contextual_control Rat_Poverty Rat_Mig_Cum Pop_Den
{txt}
{com}. 
. local religion Rat_GC_ProE Rat_GC_ProM Rat_GC_ProB Rat_GC_Jew Rat_GC_Oth
{txt}
{com}. local contextual_control Rat_Poverty Rat_Mig_Cum Pop_Den
{txt}
{com}. 
. local demographics_raw i.Female c.RAT_Female i.AgeGrp4 c.RAT_AgeGrp4_2 c.RAT_AgeGrp4_3 c.RAT_AgeGrp4_4 i.Race5 c.RAT_Race5_2 c.RAT_Race5_3 c.RAT_Race5_4 c.RAT_Race5_5 i.BornUSA c.RAT_BornUSA i.MarStat5 c.RAT_MarStat5_2 c.RAT_MarStat5_3 c.RAT_MarStat5_4 c.RAT_MarStat5_5
{txt}
{com}. local demographics_same i.Female c.std_same_prop_Sex i.AgeGrp4 c.std_same_prop_AgeGrp4 i.Race5 c.std_same_prop_Race5 i.BornUSA c.std_same_prop_BornUSA i.MarStat5 c.std_same_prop_MarStat5 
{txt}
{com}. local demographics_inter i.Female##c.std_same_prop_Sex i.AgeGrp4##c.std_same_prop_AgeGrp4 i.Race5##c.std_same_prop_Race5 i.BornUSA##c.std_same_prop_BornUSA i.MarStat5##c.std_same_prop_MarStat5
{txt}
{com}. 
. 
. if (`model_version' == 1){c -(}
.         local model_eq i.Year `demographics_raw' `contextual_control' `religion' 
.         local margin_demographics Female RAT_Female AgeGrp4 RAT_AgeGrp4_2 RAT_AgeGrp4_3 RAT_AgeGrp4_4 Race5 RAT_Race5_2 RAT_Race5_3 RAT_Race5_4 RAT_Race5_5 BornUSA RAT_BornUSA MarStat5 RAT_MarStat5_2 RAT_MarStat5_3 RAT_MarStat5_4 RAT_MarStat5_5 
. {c )-}
{txt}
{com}. if (`model_version' == 2){c -(}
.         local model_eq i.Year `demographics_raw' `contextual_control' `religion' i.UnEmpl c.RAT_UnEmpl i.PhysProb c.RAT_PhysProb
.         local margin_demographics Female RAT_Female AgeGrp4 RAT_AgeGrp4_2 RAT_AgeGrp4_3 RAT_AgeGrp4_4 Race5 RAT_Race5_2 RAT_Race5_3 RAT_Race5_4 RAT_Race5_5 BornUSA RAT_BornUSA MarStat5 RAT_MarStat5_2 RAT_MarStat5_3 RAT_MarStat5_4 RAT_MarStat5_5 UnEmpl RAT_UnEmpl PhysProb RAT_PhysProb
. {c )-}
{txt}
{com}. if (`model_version' == 3){c -(}
.         local model_eq i.Year `demographics_same' `contextual_control' `religion' 
.         local margin_demographics Female std_same_prop_Sex AgeGrp4 std_same_prop_AgeGrp4 Race5 std_same_prop_Race5 BornUSA std_same_prop_BornUSA MarStat5 std_same_prop_MarStat5 
. {c )-}
{txt}
{com}. if (`model_version' == 4){c -(}
.         local model_eq i.Year `demographics_same' `contextual_control' `religion' i.UnEmpl c.std_same_prop_UnEmpl i.PhysProb c.std_same_prop_PhysProb
.         local margin_demographics Female std_same_prop_Sex AgeGrp4 std_same_prop_AgeGrp4 Race5 std_same_prop_Race5 BornUSA std_same_prop_BornUSA MarStat5 std_same_prop_MarStat5 UnEmpl std_same_prop_UnEmpl PhysProb std_same_prop_PhysProb
. {c )-}
{txt}
{com}. if (`model_version' == 5){c -(}
.         local model_eq i.Year `demographics_inter' `contextual_control' `religion' 
. {c )-}
{txt}
{com}. if (`model_version' == 6){c -(}
.         local model_eq i.Year `demographics_inter' `contextual_control' `religion' i.UnEmpl##c.std_same_prop_UnEmpl i.PhysProb##c.std_same_prop_PhysProb
. {c )-}       
{txt}
{com}. * test how long it would take.
. * mi estimate: svy: mean Suic 
. * local demographics i.Female c.RAT_Female i.AgeGrp4 c.RAT_AgeGrp4_2 c.RAT_AgeGrp4_3 c.RAT_AgeGrp4_4 i.Race5 c.RAT_Race5_2 c.RAT_Race5_3 c.RAT_Race5_4 c.RAT_Race5_5 i.BornUSA c.RAT_BornUSA i.MarStat5 c.RAT_MarStat5_2 c.RAT_MarStat5_3 c.RAT_MarStat5_4 c.RAT_MarStat5_5
. * mi estimate: svy: logit Suic i.St `demographics' UnEmpl RAT_UnEmpl PhysProb RAT_PhysProb
. 
. * main effects : margins
. if ("`model'" == "mi"){c -(}
. 
.         mi estimate: svy: logit Suic i.St `model_eq', or 
{res}
{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Survey: Logistic regression{col 47}Number of obs{col 65}= {res}  11,814,307

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918
{txt}{col 47}Average RVI{col 65}= {res}      0.0004
{txt}{col 47}Largest FMI{col 65}= {res}      0.0043
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      909.35
{txt}{col 47}        avg{col 65}= {res}      914.60
{txt}{col 47}        max{col 65}= {res}      915.01
{txt}Model F test:{ralign 16: {res:Equal FMI}}{col 47}F({res}  55{txt},{res}  915.0{txt}){col 65}= {res}      618.69
{txt}Within VCE type: {ralign 12:{res:Linearized}}{col 47}Prob > F{col 65}= {res}      0.0000

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}        Suic{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 10}St {c |}
{space 10}8  {c |}{col 14}{res}{space 2}  -.16788{col 26}{space 2} .0971364{col 37}{space 1}   -1.73{col 46}{space 3}0.084{col 54}{space 4}-.3585161{col 67}{space 3} .0227561
{txt}{space 9}13  {c |}{col 14}{res}{space 2}-.5001434{col 26}{space 2} .1004458{col 37}{space 1}   -4.98{col 46}{space 3}0.000{col 54}{space 4}-.6972742{col 67}{space 3}-.3030125
{txt}{space 9}21  {c |}{col 14}{res}{space 2}-.5727156{col 26}{space 2} .0978785{col 37}{space 1}   -5.85{col 46}{space 3}0.000{col 54}{space 4}-.7648081{col 67}{space 3}-.3806232
{txt}{space 9}24  {c |}{col 14}{res}{space 2}-.7057739{col 26}{space 2} .1003643{col 37}{space 1}   -7.03{col 46}{space 3}0.000{col 54}{space 4}-.9027449{col 67}{space 3}-.5088029
{txt}{space 9}25  {c |}{col 14}{res}{space 2}-.8917392{col 26}{space 2} .0987201{col 37}{space 1}   -9.03{col 46}{space 3}0.000{col 54}{space 4}-1.085483{col 67}{space 3}-.6979951
{txt}{space 9}34  {c |}{col 14}{res}{space 2}-.8595632{col 26}{space 2} .1042853{col 37}{space 1}   -8.24{col 46}{space 3}0.000{col 54}{space 4}-1.064229{col 67}{space 3}-.6548971
{txt}{space 9}35  {c |}{col 14}{res}{space 2}-.2604716{col 26}{space 2}  .145413{col 37}{space 1}   -1.79{col 46}{space 3}0.074{col 54}{space 4}-.5458533{col 67}{space 3} .0249101
{txt}{space 9}37  {c |}{col 14}{res}{space 2}-.4768946{col 26}{space 2} .1021864{col 37}{space 1}   -4.67{col 46}{space 3}0.000{col 54}{space 4}-.6774415{col 67}{space 3}-.2763477
{txt}{space 9}40  {c |}{col 14}{res}{space 2} -.691179{col 26}{space 2} .1068607{col 37}{space 1}   -6.47{col 46}{space 3}0.000{col 54}{space 4}-.9008996{col 67}{space 3}-.4814584
{txt}{space 9}41  {c |}{col 14}{res}{space 2}-.3070902{col 26}{space 2} .0949924{col 37}{space 1}   -3.23{col 46}{space 3}0.001{col 54}{space 4}-.4935185{col 67}{space 3}-.1206619
{txt}{space 9}44  {c |}{col 14}{res}{space 2}-.7736412{col 26}{space 2} .0978637{col 37}{space 1}   -7.91{col 46}{space 3}0.000{col 54}{space 4}-.9657046{col 67}{space 3}-.5815779
{txt}{space 9}45  {c |}{col 14}{res}{space 2}-.4922061{col 26}{space 2} .1013778{col 37}{space 1}   -4.86{col 46}{space 3}0.000{col 54}{space 4}-.6911661{col 67}{space 3}-.2932461
{txt}{space 9}49  {c |}{col 14}{res}{space 2} .4905819{col 26}{space 2} .1663927{col 37}{space 1}    2.95{col 46}{space 3}0.003{col 54}{space 4} .1640262{col 67}{space 3} .8171376
{txt}{space 9}51  {c |}{col 14}{res}{space 2}-.3954205{col 26}{space 2} .0978763{col 37}{space 1}   -4.04{col 46}{space 3}0.000{col 54}{space 4}-.5875086{col 67}{space 3}-.2033324
{txt}{space 9}55  {c |}{col 14}{res}{space 2}-.6264339{col 26}{space 2} .0981977{col 37}{space 1}   -6.38{col 46}{space 3}0.000{col 54}{space 4}-.8191527{col 67}{space 3}-.4337151
{txt}{space 12} {c |}
{space 8}Year {c |}
{space 7}2006  {c |}{col 14}{res}{space 2}-.0530007{col 26}{space 2} .0186056{col 37}{space 1}   -2.85{col 46}{space 3}0.004{col 54}{space 4}-.0895154{col 67}{space 3} -.016486
{txt}{space 7}2007  {c |}{col 14}{res}{space 2} .0236657{col 26}{space 2} .0197091{col 37}{space 1}    1.20{col 46}{space 3}0.230{col 54}{space 4}-.0150147{col 67}{space 3}  .062346
{txt}{space 7}2008  {c |}{col 14}{res}{space 2} .0254669{col 26}{space 2} .0219947{col 37}{space 1}    1.16{col 46}{space 3}0.247{col 54}{space 4} -.017699{col 67}{space 3} .0686328
{txt}{space 7}2009  {c |}{col 14}{res}{space 2} .0799392{col 26}{space 2} .0218966{col 37}{space 1}    3.65{col 46}{space 3}0.000{col 54}{space 4} .0369658{col 67}{space 3} .1229126
{txt}{space 7}2010  {c |}{col 14}{res}{space 2} .0692682{col 26}{space 2} .0249575{col 37}{space 1}    2.78{col 46}{space 3}0.006{col 54}{space 4} .0202876{col 67}{space 3} .1182487
{txt}{space 7}2011  {c |}{col 14}{res}{space 2} .1017935{col 26}{space 2} .0276315{col 37}{space 1}    3.68{col 46}{space 3}0.000{col 54}{space 4}  .047565{col 67}{space 3}  .156022
{txt}{space 12} {c |}
{space 4}1.Female {c |}{col 14}{res}{space 2}-1.397806{col 26}{space 2} .0150897{col 37}{space 1}  -92.63{col 46}{space 3}0.000{col 54}{space 4} -1.42742{col 67}{space 3}-1.368192
{txt}{space 2}RAT_Female {c |}{col 14}{res}{space 2} .0015534{col 26}{space 2} .0055304{col 37}{space 1}    0.28{col 46}{space 3}0.779{col 54}{space 4}-.0093002{col 67}{space 3} .0124071
{txt}{space 12} {c |}
{space 5}AgeGrp4 {c |}
{space 10}2  {c |}{col 14}{res}{space 2} .7578877{col 26}{space 2} .0224473{col 37}{space 1}   33.76{col 46}{space 3}0.000{col 54}{space 4} .7138336{col 67}{space 3} .8019418
{txt}{space 10}3  {c |}{col 14}{res}{space 2} .8956439{col 26}{space 2} .0244269{col 37}{space 1}   36.67{col 46}{space 3}0.000{col 54}{space 4} .8477045{col 67}{space 3} .9435832
{txt}{space 10}4  {c |}{col 14}{res}{space 2} .6852022{col 26}{space 2} .0270129{col 37}{space 1}   25.37{col 46}{space 3}0.000{col 54}{space 4} .6321875{col 67}{space 3} .7382168
{txt}{space 12} {c |}
RAT_AgeGrp~2 {c |}{col 14}{res}{space 2}-.0062893{col 26}{space 2} .0030765{col 37}{space 1}   -2.04{col 46}{space 3}0.041{col 54}{space 4}-.0123271{col 67}{space 3}-.0002515
{txt}RAT_AgeGrp~3 {c |}{col 14}{res}{space 2}-.0153696{col 26}{space 2} .0037043{col 37}{space 1}   -4.15{col 46}{space 3}0.000{col 54}{space 4}-.0226395{col 67}{space 3}-.0080997
{txt}RAT_AgeGrp~4 {c |}{col 14}{res}{space 2} .0115652{col 26}{space 2} .0044074{col 37}{space 1}    2.62{col 46}{space 3}0.009{col 54}{space 4} .0029155{col 67}{space 3}  .020215
{txt}{space 12} {c |}
{space 7}Race5 {c |}
{space 10}2  {c |}{col 14}{res}{space 2} -1.01711{col 26}{space 2} .0260642{col 37}{space 1}  -39.02{col 46}{space 3}0.000{col 54}{space 4}-1.068263{col 67}{space 3}-.9659577
{txt}{space 10}3  {c |}{col 14}{res}{space 2} .0157149{col 26}{space 2} .0775204{col 37}{space 1}    0.20{col 46}{space 3}0.839{col 54}{space 4}-.1364235{col 67}{space 3} .1678532
{txt}{space 10}4  {c |}{col 14}{res}{space 2}-.4143794{col 26}{space 2} .0460568{col 37}{space 1}   -9.00{col 46}{space 3}0.000{col 54}{space 4}-.5047686{col 67}{space 3}-.3239901
{txt}{space 10}5  {c |}{col 14}{res}{space 2}-.7125329{col 26}{space 2} .0406941{col 37}{space 1}  -17.51{col 46}{space 3}0.000{col 54}{space 4}-.7923974{col 67}{space 3}-.6326684
{txt}{space 12} {c |}
{space 1}RAT_Race5_2 {c |}{col 14}{res}{space 2} .0041289{col 26}{space 2} .0009892{col 37}{space 1}    4.17{col 46}{space 3}0.000{col 54}{space 4} .0021875{col 67}{space 3} .0060703
{txt}{space 1}RAT_Race5_3 {c |}{col 14}{res}{space 2} .0076817{col 26}{space 2} .0034166{col 37}{space 1}    2.25{col 46}{space 3}0.025{col 54}{space 4} .0009764{col 67}{space 3}  .014387
{txt}{space 1}RAT_Race5_4 {c |}{col 14}{res}{space 2} .0059052{col 26}{space 2}  .005105{col 37}{space 1}    1.16{col 46}{space 3}0.248{col 54}{space 4}-.0041137{col 67}{space 3}  .015924
{txt}{space 1}RAT_Race5_5 {c |}{col 14}{res}{space 2}  .008181{col 26}{space 2} .0032489{col 37}{space 1}    2.52{col 46}{space 3}0.012{col 54}{space 4}  .001805{col 67}{space 3} .0145571
{txt}{space 3}1.BornUSA {c |}{col 14}{res}{space 2} .4633676{col 26}{space 2}   .03477{col 37}{space 1}   13.33{col 46}{space 3}0.000{col 54}{space 4} .3951294{col 67}{space 3} .5316058
{txt}{space 1}RAT_BornUSA {c |}{col 14}{res}{space 2} .0090352{col 26}{space 2} .0034721{col 37}{space 1}    2.60{col 46}{space 3}0.009{col 54}{space 4} .0022209{col 67}{space 3} .0158495
{txt}{space 12} {c |}
{space 4}MarStat5 {c |}
{space 10}2  {c |}{col 14}{res}{space 2} .8264811{col 26}{space 2}  .021852{col 37}{space 1}   37.82{col 46}{space 3}0.000{col 54}{space 4} .7835949{col 67}{space 3} .8693673
{txt}{space 10}3  {c |}{col 14}{res}{space 2} 1.128002{col 26}{space 2} .0149022{col 37}{space 1}   75.69{col 46}{space 3}0.000{col 54}{space 4} 1.098755{col 67}{space 3} 1.157248
{txt}{space 10}4  {c |}{col 14}{res}{space 2}-.0815195{col 26}{space 2} .0910885{col 37}{space 1}   -0.89{col 46}{space 3}0.371{col 54}{space 4} -.260287{col 67}{space 3}  .097248
{txt}{space 10}5  {c |}{col 14}{res}{space 2} .8075421{col 26}{space 2} .0178729{col 37}{space 1}   45.18{col 46}{space 3}0.000{col 54}{space 4} .7724653{col 67}{space 3}  .842619
{txt}{space 12} {c |}
RAT_MarSta~2 {c |}{col 14}{res}{space 2}-.0099811{col 26}{space 2} .0074974{col 37}{space 1}   -1.33{col 46}{space 3}0.183{col 54}{space 4}-.0246953{col 67}{space 3} .0047331
{txt}RAT_MarSta~3 {c |}{col 14}{res}{space 2} .0058123{col 26}{space 2} .0044998{col 37}{space 1}    1.29{col 46}{space 3}0.197{col 54}{space 4}-.0030189{col 67}{space 3} .0146435
{txt}RAT_MarSta~4 {c |}{col 14}{res}{space 2} .0057564{col 26}{space 2} .0090382{col 37}{space 1}    0.64{col 46}{space 3}0.524{col 54}{space 4}-.0119817{col 67}{space 3} .0234945
{txt}RAT_MarSta~5 {c |}{col 14}{res}{space 2}-.0084171{col 26}{space 2} .0024424{col 37}{space 1}   -3.45{col 46}{space 3}0.001{col 54}{space 4}-.0132104{col 67}{space 3}-.0036238
{txt}{space 1}Rat_Poverty {c |}{col 14}{res}{space 2}-.0023226{col 26}{space 2} .0020953{col 37}{space 1}   -1.11{col 46}{space 3}0.268{col 54}{space 4}-.0064347{col 67}{space 3} .0017895
{txt}{space 1}Rat_Mig_Cum {c |}{col 14}{res}{space 2}-1.630236{col 26}{space 2} .4607926{col 37}{space 1}   -3.54{col 46}{space 3}0.000{col 54}{space 4}-2.534569{col 67}{space 3} -.725903
{txt}{space 5}Pop_Den {c |}{col 14}{res}{space 2} 4.78e-06{col 26}{space 2} 4.94e-06{col 37}{space 1}    0.97{col 46}{space 3}0.333{col 54}{space 4}-4.91e-06{col 67}{space 3} .0000145
{txt}{space 1}Rat_GC_ProE {c |}{col 14}{res}{space 2} 5.16e-06{col 26}{space 2} 8.92e-06{col 37}{space 1}    0.58{col 46}{space 3}0.563{col 54}{space 4}-.0000124{col 67}{space 3} .0000227
{txt}{space 1}Rat_GC_ProM {c |}{col 14}{res}{space 2} 4.93e-06{col 26}{space 2} .0000112{col 37}{space 1}    0.44{col 46}{space 3}0.659{col 54}{space 4} -.000017{col 67}{space 3} .0000268
{txt}{space 1}Rat_GC_ProB {c |}{col 14}{res}{space 2}-.0000593{col 26}{space 2} .0000421{col 37}{space 1}   -1.41{col 46}{space 3}0.160{col 54}{space 4}-.0001419{col 67}{space 3} .0000234
{txt}{space 2}Rat_GC_Jew {c |}{col 14}{res}{space 2} .0000893{col 26}{space 2}  .000075{col 37}{space 1}    1.19{col 46}{space 3}0.234{col 54}{space 4}-.0000578{col 67}{space 3} .0002365
{txt}{space 2}Rat_GC_Oth {c |}{col 14}{res}{space 2}-.0001056{col 26}{space 2} .0000229{col 37}{space 1}   -4.62{col 46}{space 3}0.000{col 54}{space 4}-.0001505{col 67}{space 3}-.0000608
{txt}{space 7}_cons {c |}{col 14}{res}{space 2} -9.35924{col 26}{space 2} .5490225{col 37}{space 1}  -17.05{col 46}{space 3}0.000{col 54}{space 4}-10.43673{col 67}{space 3} -8.28175
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}{com}.         estimates store m1 
. 
.         if (`model_version' <= 4){c -(}
.                 estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                 mchange_mi Female "binary" "mi"

                           {txt}Female
{hline 61}
      Percentiles      Smallest
 1%    {res}        0              0
{txt} 5%    {res}        0              0
{txt}10%    {res}        0              0       {txt}Obs         {res}135,582,733
{txt}25%    {res}        0              0       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res}        1                      {txt}Mean          {res} .5158124
                        {txt}Largest       Std. Dev.     {res} .4997499
{txt}75%    {res}        1              1
{txt}90%    {res}        1              1       {txt}Variance      {res}   .24975
{txt}95%    {res}        1              1       {txt}Skewness      {res}-.0632814
{txt}99%    {res}        1              1       {txt}Kurtosis      {res} 1.004005

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0001
{txt}{col 47}Largest FMI{col 65}= {res}      0.0001
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      914.92
{txt}{col 47}        avg{col 65}= {res}      914.95
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      914.98

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:Female}{space 10}{txt:=} {space 10}0}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:Female}{space 10}{txt:=} {space 10}1}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0002464{col 26}{space 2} 1.85e-06{col 37}{space 1}  133.31{col 46}{space 3}0.000{col 54}{space 4} .0002428{col 67}{space 3} .0002501
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0000609{col 26}{space 2} 7.93e-07{col 37}{space 1}   76.79{col 46}{space 3}0.000{col 54}{space 4} .0000594{col 67}{space 3} .0000625
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:-0.0001855}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000020}}}{space 1}{space 1}{ralign 8:{res:{sf:-9.2e+01}}}{space 1}{space 1}{ralign 8:{res:{sf: 0.0e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0001895}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0001816}}}{space 1}
{com}.                 estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                 mchange_mi AgeGrp4 "categorical" "mi"
{res}
{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Pairwise comparisons of predictive margins{col 47}Number of obs{col 65}= {res}  11,079,265

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0006
{txt}{col 47}Largest FMI{col 65}= {res}      0.0012
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      913.80
{txt}{col 47}        avg{col 65}= {res}      914.25
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      914.96

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 14}{hline 12}
{col 14}{c |}   Contrast{col 26}   Std. Err.{col 38}     [95% Con{col 51}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 14}{hline 12}
{space 5}AgeGrp4 {c |}
{space 5}2 vs 1  {c |}{col 14}{res}{space 2}  .000087{col 26}{space 2} 2.30e-06{col 37}{space 5} .0000824{col 51}{space 3} .0000915
{txt}{space 5}3 vs 1  {c |}{col 14}{res}{space 2} .0001111{col 26}{space 2} 2.68e-06{col 37}{space 5} .0001059{col 51}{space 3} .0001164
{txt}{space 5}4 vs 1  {c |}{col 14}{res}{space 2} .0000755{col 26}{space 2} 2.93e-06{col 37}{space 5} .0000697{col 51}{space 3} .0000812
{txt}{space 5}3 vs 2  {c |}{col 14}{res}{space 2} .0000242{col 26}{space 2} 2.44e-06{col 37}{space 5} .0000194{col 51}{space 3}  .000029
{txt}{space 5}4 vs 2  {c |}{col 14}{res}{space 2}-.0000115{col 26}{space 2} 3.05e-06{col 37}{space 5}-.0000174{col 51}{space 3}-5.49e-06
{txt}{space 5}4 vs 3  {c |}{col 14}{res}{space 2}-.0000356{col 26}{space 2} 2.83e-06{col 37}{space 5}-.0000412{col 51}{space 3}-.0000301
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 14}{hline 12}
{res}{com}.                 estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                 mchange_mi Race5 "categorical" "mi"
{res}
{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Pairwise comparisons of predictive margins{col 47}Number of obs{col 65}= {res}  11,079,265

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0000
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      914.97
{txt}{col 47}        avg{col 65}= {res}      914.99
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      915.01

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 14}{hline 12}
{col 14}{c |}   Contrast{col 26}   Std. Err.{col 38}     [95% Con{col 51}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 14}{hline 12}
{space 7}Race5 {c |}
{space 5}2 vs 1  {c |}{col 14}{res}{space 2}-.0001108{col 26}{space 2} 2.09e-06{col 37}{space 5}-.0001149{col 51}{space 3}-.0001066
{txt}{space 5}3 vs 1  {c |}{col 14}{res}{space 2} 2.75e-06{col 26}{space 2} .0000137{col 37}{space 5}-.0000241{col 51}{space 3} .0000296
{txt}{space 5}4 vs 1  {c |}{col 14}{res}{space 2}-.0000589{col 26}{space 2} 5.38e-06{col 37}{space 5}-.0000694{col 51}{space 3}-.0000483
{txt}{space 5}5 vs 1  {c |}{col 14}{res}{space 2}-.0000884{col 26}{space 2} 3.71e-06{col 37}{space 5}-.0000957{col 51}{space 3}-.0000811
{txt}{space 5}3 vs 2  {c |}{col 14}{res}{space 2} .0001135{col 26}{space 2} .0000137{col 37}{space 5} .0000867{col 51}{space 3} .0001403
{txt}{space 5}4 vs 2  {c |}{col 14}{res}{space 2} .0000519{col 26}{space 2} 5.32e-06{col 37}{space 5} .0000415{col 51}{space 3} .0000623
{txt}{space 5}5 vs 2  {c |}{col 14}{res}{space 2} .0000223{col 26}{space 2} 3.61e-06{col 37}{space 5} .0000153{col 51}{space 3} .0000294
{txt}{space 5}4 vs 3  {c |}{col 14}{res}{space 2}-.0000616{col 26}{space 2} .0000148{col 37}{space 5}-.0000906{col 51}{space 3}-.0000326
{txt}{space 5}5 vs 3  {c |}{col 14}{res}{space 2}-.0000912{col 26}{space 2} .0000145{col 37}{space 5}-.0001195{col 51}{space 3}-.0000628
{txt}{space 5}5 vs 4  {c |}{col 14}{res}{space 2}-.0000296{col 26}{space 2} 5.81e-06{col 37}{space 5} -.000041{col 51}{space 3}-.0000181
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 14}{hline 12}
{res}{com}.                 estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                 mchange_mi BornUSA "binary" "mi"

                           {txt}BornUSA
{hline 61}
      Percentiles      Smallest
 1%    {res}        0              0
{txt} 5%    {res}        0              0
{txt}10%    {res}        1              0       {txt}Obs         {res}135,582,733
{txt}25%    {res}        1              0       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res}        1                      {txt}Mean          {res} .9265123
                        {txt}Largest       Std. Dev.     {res} .2609353
{txt}75%    {res}        1              1
{txt}90%    {res}        1              1       {txt}Variance      {res} .0680872
{txt}95%    {res}        1              1       {txt}Skewness      {res}-3.269105
{txt}99%    {res}        1              1       {txt}Kurtosis      {res} 11.68704

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0000
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      915.00
{txt}{col 47}        avg{col 65}= {res}      915.00
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      915.01

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:BornUSA}{space 9}{txt:=} {space 10}0}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:BornUSA}{space 9}{txt:=} {space 10}1}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0000965{col 26}{space 2} 3.06e-06{col 37}{space 1}   31.49{col 46}{space 3}0.000{col 54}{space 4} .0000905{col 67}{space 3} .0001025
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001533{col 26}{space 2} 1.12e-06{col 37}{space 1}  137.41{col 46}{space 3}0.000{col 54}{space 4} .0001511{col 67}{space 3} .0001555
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:0.0000568}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000036}}}{space 1}{space 1}{ralign 8:{res:{sf: 1.6e+01}}}{space 1}{space 1}{ralign 8:{res:{sf: 0.0e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000498}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000639}}}{space 1}
{com}.                 estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                 mchange_mi MarStat5 "categorical" "mi"
{res}
{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Pairwise comparisons of predictive margins{col 47}Number of obs{col 65}= {res}  11,079,265

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0032
{txt}{col 47}Largest FMI{col 65}= {res}      0.0046
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      908.81
{txt}{col 47}        avg{col 65}= {res}      911.03
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      912.96

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 14}{hline 12}
{col 14}{c |}   Contrast{col 26}   Std. Err.{col 38}     [95% Con{col 51}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 14}{hline 12}
{space 4}MarStat5 {c |}
{space 5}2 vs 1  {c |}{col 14}{res}{space 2}  .000121{col 26}{space 2} 4.37e-06{col 37}{space 5} .0001124{col 51}{space 3} .0001295
{txt}{space 5}3 vs 1  {c |}{col 14}{res}{space 2} .0001966{col 26}{space 2} 3.39e-06{col 37}{space 5}   .00019{col 51}{space 3} .0002033
{txt}{space 5}4 vs 1  {c |}{col 14}{res}{space 2}-7.37e-06{col 26}{space 2} 7.93e-06{col 37}{space 5}-.0000229{col 51}{space 3} 8.19e-06
{txt}{space 5}5 vs 1  {c |}{col 14}{res}{space 2} .0001169{col 26}{space 2} 2.88e-06{col 37}{space 5} .0001113{col 51}{space 3} .0001226
{txt}{space 5}3 vs 2  {c |}{col 14}{res}{space 2} .0000757{col 26}{space 2} 5.09e-06{col 37}{space 5} .0000657{col 51}{space 3} .0000857
{txt}{space 5}4 vs 2  {c |}{col 14}{res}{space 2}-.0001283{col 26}{space 2} 8.88e-06{col 37}{space 5}-.0001458{col 51}{space 3}-.0001109
{txt}{space 5}5 vs 2  {c |}{col 14}{res}{space 2}-4.03e-06{col 26}{space 2} 5.28e-06{col 37}{space 5}-.0000144{col 51}{space 3} 6.33e-06
{txt}{space 5}4 vs 3  {c |}{col 14}{res}{space 2} -.000204{col 26}{space 2} 7.99e-06{col 37}{space 5}-.0002197{col 51}{space 3}-.0001883
{txt}{space 5}5 vs 3  {c |}{col 14}{res}{space 2}-.0000797{col 26}{space 2} 4.16e-06{col 37}{space 5}-.0000879{col 51}{space 3}-.0000715
{txt}{space 5}5 vs 4  {c |}{col 14}{res}{space 2} .0001243{col 26}{space 2} 7.80e-06{col 37}{space 5}  .000109{col 51}{space 3} .0001396
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 14}{hline 12}
{res}{com}.                 
.                 if (`model_version' == 1 | `model_version' == 2) {c -(}
.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi RAT_Female "continuous" "mi"

                         {txt}RAT_Female
{hline 61}
      Percentiles      Smallest
 1%    {res} 45.39424       43.22342
{txt} 5%    {res} 48.12819       43.22342
{txt}10%    {res} 49.38633       43.22342       {txt}Obs         {res}135,582,733
{txt}25%    {res} 50.19465       43.22342       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res} 51.19414                      {txt}Mean          {res} 51.04837
                        {txt}Largest       Std. Dev.     {res} 1.607091
{txt}75%    {res} 52.14429       56.04515
{txt}90%    {res} 52.77465       56.04515       {txt}Variance      {res} 2.582741
{txt}95%    {res} 53.17716       56.04515       {txt}Skewness      {res}-1.062898
{txt}99%    {res}  54.2445       56.04515       {txt}Kurtosis      {res} 5.535048

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0000
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      915.01
{txt}{col 47}        avg{col 65}= {res}      915.01
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      915.01

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:RAT_Female}{space 6}{txt:=} {space 3}43.22342}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:RAT_Female}{space 6}{txt:=} {space 3}56.04515}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001466{col 26}{space 2} 6.56e-06{col 37}{space 1}   22.36{col 46}{space 3}0.000{col 54}{space 4} .0001337{col 67}{space 3} .0001595
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001496{col 26}{space 2} 4.10e-06{col 37}{space 1}   36.48{col 46}{space 3}0.000{col 54}{space 4} .0001415{col 67}{space 3} .0001576
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:0.0000029}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000105}}}{space 1}{space 1}{ralign 8:{res:{sf:0.2815796}}}{space 1}{space 1}{ralign 8:{res:{sf:0.7783294}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000176}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000235}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi RAT_AgeGrp4_2 "continuous" "mi"

                        {txt}RAT_AgeGrp4_2
{hline 61}
      Percentiles      Smallest
 1%    {res} 24.41741       18.69043
{txt} 5%    {res} 26.25468       18.69043
{txt}10%    {res} 27.48339       18.69043       {txt}Obs         {res}135,582,733
{txt}25%    {res} 29.49091       18.69043       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res} 32.05476                      {txt}Mean          {res} 32.44892
                        {txt}Largest       Std. Dev.     {res} 4.216616
{txt}75%    {res} 34.66521       51.79797
{txt}90%    {res} 38.41551       51.79797       {txt}Variance      {res} 17.77985
{txt}95%    {res} 40.72519       51.79797       {txt}Skewness      {res} .6572291
{txt}99%    {res} 44.04443       51.79797       {txt}Kurtosis      {res} 3.661787

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0000
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      915.00
{txt}{col 47}        avg{col 65}= {res}      915.00
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      915.00

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:RAT_AgeGrp~2}{space 4}{txt:=} {space 3}18.69043}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:RAT_AgeGrp~2}{space 4}{txt:=} {space 3}51.79797}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001632{col 26}{space 2} 7.66e-06{col 37}{space 1}   21.31{col 46}{space 3}0.000{col 54}{space 4} .0001482{col 67}{space 3} .0001782
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001325{col 26}{space 2} 7.39e-06{col 37}{space 1}   17.94{col 46}{space 3}0.000{col 54}{space 4}  .000118{col 67}{space 3}  .000147
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:-0.0000307}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000149}}}{space 1}{space 1}{ralign 8:{res:{sf:-2.1e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0401255}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000600}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000014}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi RAT_AgeGrp4_3 "continuous" "mi"

                        {txt}RAT_AgeGrp4_3
{hline 61}
      Percentiles      Smallest
 1%    {res} 23.56008       19.29798
{txt} 5%    {res} 27.30353       19.29798
{txt}10%    {res} 28.88125       19.29798       {txt}Obs         {res}135,582,733
{txt}25%    {res}  30.6143       19.29798       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res} 32.15104                      {txt}Mean          {res} 32.09528
                        {txt}Largest       Std. Dev.     {res} 2.802824
{txt}75%    {res}   33.866        47.4186
{txt}90%    {res} 35.46149        47.4186       {txt}Variance      {res} 7.855821
{txt}95%    {res}  36.2629        47.4186       {txt}Skewness      {res}-.4997261
{txt}99%    {res}  38.0363        47.4186       {txt}Kurtosis      {res} 4.789628

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0000
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      915.00
{txt}{col 47}        avg{col 65}= {res}      915.00
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      915.00

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:RAT_AgeGrp~3}{space 4}{txt:=} {space 3}19.29798}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:RAT_AgeGrp~3}{space 4}{txt:=} {space 4}47.4186}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001803{col 26}{space 2} 8.69e-06{col 37}{space 1}   20.74{col 46}{space 3}0.000{col 54}{space 4} .0001632{col 67}{space 3} .0001973
{txt}{space 10}2  {c |}{col 14}{res}{space 2}  .000117{col 26}{space 2} 6.64e-06{col 37}{space 1}   17.63{col 46}{space 3}0.000{col 54}{space 4}  .000104{col 67}{space 3} .0001301
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:-0.0000633}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000152}}}{space 1}{space 1}{ralign 8:{res:{sf:-4.2e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000352}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000931}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000334}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi RAT_AgeGrp4_4 "continuous" "mi"

                        {txt}RAT_AgeGrp4_4
{hline 61}
      Percentiles      Smallest
 1%    {res} 8.643331       6.363428
{txt} 5%    {res} 10.62395       6.363428
{txt}10%    {res} 12.01889       6.363428       {txt}Obs         {res}135,582,733
{txt}25%    {res} 15.03323       6.363428       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res}   17.501                      {txt}Mean          {res}   17.318
                        {txt}Largest       Std. Dev.     {res} 3.740434
{txt}75%    {res} 20.01174       29.62251
{txt}90%    {res} 21.92176       29.62251       {txt}Variance      {res} 13.99085
{txt}95%    {res} 22.92662       29.62251       {txt}Skewness      {res}-.1965066
{txt}99%    {res} 25.70682       29.62251       {txt}Kurtosis      {res} 2.852231

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0000
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      915.00
{txt}{col 47}        avg{col 65}= {res}      915.00
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      915.00

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:RAT_AgeGrp~4}{space 4}{txt:=} {space 3}6.363428}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:RAT_AgeGrp~4}{space 4}{txt:=} {space 3}29.62251}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001326{col 26}{space 2} 5.78e-06{col 37}{space 1}   22.93{col 46}{space 3}0.000{col 54}{space 4} .0001212{col 67}{space 3} .0001439
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001735{col 26}{space 2} .0000104{col 37}{space 1}   16.75{col 46}{space 3}0.000{col 54}{space 4} .0001532{col 67}{space 3} .0001938
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:0.0000409}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000160}}}{space 1}{space 1}{ralign 8:{res:{sf: 2.6e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0108115}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000095}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000723}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi RAT_Race5_2 "continuous" "mi"

                         {txt}RAT_Race5_2
{hline 61}
      Percentiles      Smallest
 1%    {res} .0990109              0
{txt} 5%    {res} .3381911              0
{txt}10%    {res} .6039064              0       {txt}Obs         {res}135,582,733
{txt}25%    {res}  1.59464              0       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res}   5.7945                      {txt}Mean          {res} 12.30798
                        {txt}Largest       Std. Dev.     {res} 14.12576
{txt}75%    {res} 20.08249       65.45666
{txt}90%    {res} 33.72202       65.45666       {txt}Variance      {res}  199.537
{txt}95%    {res} 42.64764       65.45666       {txt}Skewness      {res} 1.379799
{txt}99%    {res} 56.82112       65.45666       {txt}Kurtosis      {res} 4.208692

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0000
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      915.00
{txt}{col 47}        avg{col 65}= {res}      915.00
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      915.00

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:RAT_Race5_2}{space 5}{txt:=} {space 10}0}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:RAT_Race5_2}{space 5}{txt:=} {space 3}65.45666}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001409{col 26}{space 2} 1.81e-06{col 37}{space 1}   77.98{col 46}{space 3}0.000{col 54}{space 4} .0001374{col 67}{space 3} .0001445
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001846{col 26}{space 2} 9.97e-06{col 37}{space 1}   18.52{col 46}{space 3}0.000{col 54}{space 4} .0001651{col 67}{space 3} .0002042
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:0.0000437}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000115}}}{space 1}{space 1}{ralign 8:{res:{sf: 3.8e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0001474}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000212}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000662}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi RAT_Race5_3 "continuous" "mi"

                         {txt}RAT_Race5_3
{hline 61}
      Percentiles      Smallest
 1%    {res}        0              0
{txt} 5%    {res} .0018458              0
{txt}10%    {res} .0539814              0       {txt}Obs         {res}135,582,733
{txt}25%    {res} .1379112              0       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res} .3117765                      {txt}Mean          {res} 1.842507
                        {txt}Largest       Std. Dev.     {res} 5.832106
{txt}75%    {res} 1.009428       64.16144
{txt}90%    {res} 4.032134       64.16144       {txt}Variance      {res} 34.01346
{txt}95%    {res} 7.956508       64.16144       {txt}Skewness      {res} 7.046266
{txt}99%    {res} 48.06156       64.16144       {txt}Kurtosis      {res} 58.98451

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0000
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      915.01
{txt}{col 47}        avg{col 65}= {res}      915.01
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      915.01

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:RAT_Race5_3}{space 5}{txt:=} {space 10}0}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:RAT_Race5_3}{space 5}{txt:=} {space 3}64.16144}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2}  .000147{col 26}{space 2} 1.09e-06{col 37}{space 1}  134.76{col 46}{space 3}0.000{col 54}{space 4} .0001449{col 67}{space 3} .0001491
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0002406{col 26}{space 2} .0000519{col 37}{space 1}    4.64{col 46}{space 3}0.000{col 54}{space 4} .0001388{col 67}{space 3} .0003424
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:0.0000936}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000524}}}{space 1}{space 1}{ralign 8:{res:{sf: 1.8e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0742602}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000092}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0001964}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi RAT_Race5_4 "continuous" "mi"

                         {txt}RAT_Race5_4
{hline 61}
      Percentiles      Smallest
 1%    {res}        0              0
{txt} 5%    {res} .1463881              0
{txt}10%    {res} .2372605              0       {txt}Obs         {res}135,582,733
{txt}25%    {res} .4961168              0       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res} 1.117422                      {txt}Mean          {res} 2.200708
                        {txt}Largest       Std. Dev.     {res} 2.988983
{txt}75%    {res}  2.52559       21.00913
{txt}90%    {res}  5.15511       21.00913       {txt}Variance      {res} 8.934017
{txt}95%    {res}  8.65602       21.00913       {txt}Skewness      {res} 2.837257
{txt}99%    {res} 15.57015       21.00913       {txt}Kurtosis      {res}  12.5286

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0000
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      915.01
{txt}{col 47}        avg{col 65}= {res}      915.01
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      915.01

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:RAT_Race5_4}{space 5}{txt:=} {space 10}0}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:RAT_Race5_4}{space 5}{txt:=} {space 3}21.00913}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001459{col 26}{space 2} 2.29e-06{col 37}{space 1}   63.80{col 46}{space 3}0.000{col 54}{space 4} .0001414{col 67}{space 3} .0001504
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001651{col 26}{space 2} .0000154{col 37}{space 1}   10.73{col 46}{space 3}0.000{col 54}{space 4}  .000135{col 67}{space 3} .0001953
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:0.0000193}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000174}}}{space 1}{space 1}{ralign 8:{res:{sf: 1.1e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.2694268}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000149}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000535}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi RAT_Race5_5 "continuous" "mi"

                         {txt}RAT_Race5_5
{hline 61}
      Percentiles      Smallest
 1%    {res} .3568469              0
{txt} 5%    {res}  .795982              0
{txt}10%    {res}  1.14643              0       {txt}Obs         {res}135,582,733
{txt}25%    {res} 2.027097              0       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res} 4.047258                      {txt}Mean          {res} 6.821887
                        {txt}Largest       Std. Dev.     {res}  8.43017
{txt}75%    {res} 8.110949       71.36269
{txt}90%    {res} 14.99624       71.36269       {txt}Variance      {res} 71.06776
{txt}95%    {res} 19.11328       71.36269       {txt}Skewness      {res} 3.211222
{txt}99%    {res}   41.891       71.36269       {txt}Kurtosis      {res} 16.48395

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0000
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      915.01
{txt}{col 47}        avg{col 65}= {res}      915.01
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      915.01

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:RAT_Race5_5}{space 5}{txt:=} {space 10}0}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:RAT_Race5_5}{space 5}{txt:=} {space 3}71.36269}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001387{col 26}{space 2} 3.62e-06{col 37}{space 1}   38.28{col 46}{space 3}0.000{col 54}{space 4} .0001316{col 67}{space 3} .0001459
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0002487{col 26}{space 2} .0000514{col 37}{space 1}    4.84{col 46}{space 3}0.000{col 54}{space 4} .0001479{col 67}{space 3} .0003495
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:0.0001099}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000548}}}{space 1}{space 1}{ralign 8:{res:{sf: 2.0e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0453089}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000023}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0002176}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi RAT_BornUSA "continuous" "mi"

                         {txt}RAT_BornUSA
{hline 61}
      Percentiles      Smallest
 1%    {res} 62.70366       47.06422
{txt} 5%    {res} 75.50158       47.06422
{txt}10%    {res} 81.22276       47.06422       {txt}Obs         {res}135,582,733
{txt}25%    {res}   89.276       47.06422       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res} 93.94516                      {txt}Mean          {res} 91.46602
                        {txt}Largest       Std. Dev.     {res} 7.987417
{txt}75%    {res} 96.89728        99.7309
{txt}90%    {res} 98.08357        99.7309       {txt}Variance      {res} 63.79883
{txt}95%    {res}  98.5194        99.7309       {txt}Skewness      {res}-2.036025
{txt}99%    {res} 99.26265        99.7309       {txt}Kurtosis      {res} 7.747763

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0000
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      915.01
{txt}{col 47}        avg{col 65}= {res}      915.01
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      915.01

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:RAT_BornUSA}{space 5}{txt:=} {space 3}47.06422}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:RAT_BornUSA}{space 5}{txt:=} {space 4}99.7309}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2}  .000102{col 26}{space 2} .0000145{col 37}{space 1}    7.01{col 46}{space 3}0.000{col 54}{space 4} .0000735{col 67}{space 3} .0001306
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001642{col 26}{space 2} 6.71e-06{col 37}{space 1}   24.46{col 46}{space 3}0.000{col 54}{space 4}  .000151{col 67}{space 3} .0001773
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:0.0000621}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000212}}}{space 1}{space 1}{ralign 8:{res:{sf: 2.9e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0034025}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000206}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0001037}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi RAT_MarStat5_2 "continuous" "mi"

                       {txt}RAT_MarStat5_2
{hline 61}
      Percentiles      Smallest
 1%    {res} 2.957906       1.714519
{txt} 5%    {res} 3.972756       1.714519
{txt}10%    {res} 4.509835       1.714519       {txt}Obs         {res}135,582,733
{txt}25%    {res} 5.559734       1.714519       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res} 6.761515                      {txt}Mean          {res} 6.704567
                        {txt}Largest       Std. Dev.     {res} 1.648667
{txt}75%    {res} 7.841749       11.84405
{txt}90%    {res} 8.780593       11.84405       {txt}Variance      {res} 2.718104
{txt}95%    {res} 9.424637       11.84405       {txt}Skewness      {res}-.0327279
{txt}99%    {res} 10.50469       11.84405       {txt}Kurtosis      {res}  2.77525

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0000
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      915.00
{txt}{col 47}        avg{col 65}= {res}      915.00
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      915.00

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:RAT_MarSta~2}{space 4}{txt:=} {space 3}1.714519}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:RAT_MarSta~2}{space 4}{txt:=} {space 3}11.84405}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001554{col 26}{space 2} 5.43e-06{col 37}{space 1}   28.61{col 46}{space 3}0.000{col 54}{space 4} .0001447{col 67}{space 3}  .000166
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001404{col 26}{space 2} 5.91e-06{col 37}{space 1}   23.77{col 46}{space 3}0.000{col 54}{space 4} .0001288{col 67}{space 3}  .000152
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:-0.0000149}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000112}}}{space 1}{space 1}{ralign 8:{res:{sf:-1.3e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.1817197}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000369}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000070}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi RAT_MarStat5_3 "continuous" "mi"

                       {txt}RAT_MarStat5_3
{hline 61}
      Percentiles      Smallest
 1%    {res} 6.891079       5.235612
{txt} 5%    {res} 7.722857       5.235612
{txt}10%    {res} 8.234858       5.235612       {txt}Obs         {res}135,582,733
{txt}25%    {res} 9.401562       5.235612       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res} 10.65189                      {txt}Mean          {res} 10.74006
                        {txt}Largest       Std. Dev.     {res} 1.914127
{txt}75%    {res} 12.08289       18.81139
{txt}90%    {res} 13.33138       18.81139       {txt}Variance      {res}  3.66388
{txt}95%    {res} 13.87485       18.81139       {txt}Skewness      {res} .2119803
{txt}99%    {res} 15.27534       18.81139       {txt}Kurtosis      {res} 2.810613

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0000
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      915.00
{txt}{col 47}        avg{col 65}= {res}      915.00
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      915.00

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:RAT_MarSta~3}{space 4}{txt:=} {space 3}5.235612}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:RAT_MarSta~3}{space 4}{txt:=} {space 3}18.81139}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001437{col 26}{space 2} 3.75e-06{col 37}{space 1}   38.38{col 46}{space 3}0.000{col 54}{space 4} .0001364{col 67}{space 3} .0001511
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001555{col 26}{space 2} 5.66e-06{col 37}{space 1}   27.46{col 46}{space 3}0.000{col 54}{space 4} .0001444{col 67}{space 3} .0001667
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:0.0000118}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000092}}}{space 1}{space 1}{ralign 8:{res:{sf: 1.3e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.2001030}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000063}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000299}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi RAT_MarStat5_4 "continuous" "mi"

                       {txt}RAT_MarStat5_4
{hline 61}
      Percentiles      Smallest
 1%    {res} .5794714       .2517118
{txt} 5%    {res} .8603528       .2517118
{txt}10%    {res} 1.034397       .2517118       {txt}Obs         {res}135,582,733
{txt}25%    {res} 1.508482       .2517118       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res} 2.086036                      {txt}Mean          {res} 2.256889
                        {txt}Largest       Std. Dev.     {res} 1.031968
{txt}75%    {res} 2.864548       6.998906
{txt}90%    {res}  3.65125       6.998906       {txt}Variance      {res} 1.064958
{txt}95%    {res} 4.173747       6.998906       {txt}Skewness      {res} .8185424
{txt}99%    {res} 5.279278       6.998906       {txt}Kurtosis      {res} 3.722398

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0000
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      915.00
{txt}{col 47}        avg{col 65}= {res}      915.00
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      915.00

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:RAT_MarSta~4}{space 4}{txt:=} {space 3}.2517118}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:RAT_MarSta~4}{space 4}{txt:=} {space 3}6.998906}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001468{col 26}{space 2} 2.95e-06{col 37}{space 1}   49.68{col 46}{space 3}0.000{col 54}{space 4}  .000141{col 67}{space 3} .0001526
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001526{col 26}{space 2} 6.47e-06{col 37}{space 1}   23.60{col 46}{space 3}0.000{col 54}{space 4} .0001399{col 67}{space 3} .0001653
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:0.0000058}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000092}}}{space 1}{space 1}{ralign 8:{res:{sf:0.6322142}}}{space 1}{space 1}{ralign 8:{res:{sf:0.5274045}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000122}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000238}}}{space 1}
{com}.                         estimates restore m1
{txt}(results {stata estimates replay m1:m1} are active now)
{com}.                         mchange_mi RAT_MarStat5_5 "continuous" "mi"

                       {txt}RAT_MarStat5_5
{hline 61}
      Percentiles      Smallest
 1%    {res} 19.33292       16.65138
{txt} 5%    {res} 21.06569       16.65138
{txt}10%    {res} 22.19825       16.65138       {txt}Obs         {res}135,582,733
{txt}25%    {res} 24.23431       16.65138       {txt}Sum of Wgt. {res}  135582733

{txt}50%    {res} 27.18444                      {txt}Mean          {res} 28.14833
                        {txt}Largest       Std. Dev.     {res} 5.562704
{txt}75%    {res} 30.88505        55.3223
{txt}90%    {res} 35.50837        55.3223       {txt}Variance      {res} 30.94368
{txt}95%    {res} 38.80957        55.3223       {txt}Skewness      {res} 1.157686
{txt}99%    {res} 45.98981        55.3223       {txt}Kurtosis      {res}  5.05425

{txt}Multiple-imputation estimates{col 47}Imputations{col 65}= {res}          10
{txt}Predictive margins{col 47}Number of obs{col 65}= {res}  11,764,378

{txt}{col 1}Number of strata{col 19}= {res}        1{txt}{col 47}Population size{col 65}={res}  418,786,443
{txt}{col 1}Number of PSUs{col 19}= {res}      918{txt}{col 47}Subpop. no. obs{col 65}={res}   11,079,265
{txt}{col 47}Subpop. size{col 65}={res}            .
{txt}{col 47}Average RVI{col 65}= {res}      0.0000
{txt}{col 47}Largest FMI{col 65}= {res}      0.0000
{txt}{col 47}Complete DF{col 65}= {res}         917
{txt}DF adjustment:{ralign 15: {res:Small sample}}{col 47}DF:     min{col 65}= {res}      914.98
{txt}{col 47}        avg{col 65}= {res}      914.98
{txt}Within VCE type: {ralign 12:{res:Delta-method}}{col 47}        max{col 65}= {res}      914.98

{txt}Expression{col 14}: {res}Pr(Suic), predict(pr)

{txt}{p2colset 1 14 16 2}{...}
{p2col:1._at}:{space 1}{res:{txt:RAT_MarSta~5}{space 4}{txt:=} {space 3}16.65138}{p_end}
{p2colreset}{...}

{txt}{p2colset 1 14 16 2}{...}
{p2col:2._at}:{space 1}{res:{txt:RAT_MarSta~5}{space 4}{txt:=} {space 4}55.3223}{p_end}
{p2colreset}{...}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}     Margin{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 9}_at {c |}
{space 10}1  {c |}{col 14}{res}{space 2} .0001655{col 26}{space 2} 5.48e-06{col 37}{space 1}   30.21{col 46}{space 3}0.000{col 54}{space 4} .0001548{col 67}{space 3} .0001763
{txt}{space 10}2  {c |}{col 14}{res}{space 2} .0001195{col 26}{space 2} 7.45e-06{col 37}{space 1}   16.05{col 46}{space 3}0.000{col 54}{space 4} .0001049{col 67}{space 3} .0001342
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 8:lincom}{space 1}{space 1}{ralign 8:se}{space 1}{space 1}{ralign 8:zvalue}{space 1}{space 1}{ralign 8:pvalue}{space 1}{space 1}{ralign 8:ll}{space 1}{space 1}{ralign 8:ul}{space 1}
{space 0}{hline 13}{c   +}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}{hline 10}
{space 0}{space 0}{ralign 12:1}{space 1}{c |}{space 1}{ralign 8:{res:{sf:-0.0000460}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0000128}}}{space 1}{space 1}{ralign 8:{res:{sf:-3.6e+00}}}{space 1}{space 1}{ralign 8:{res:{sf:0.0003399}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000711}}}{space 1}{space 1}{ralign 8:{res:{sf:-0.0000209}}}{space 1}
{com}.                 {c )-}
.                 if (`model_version' == 3 | `model_version' == 4) {c -(}
.                         estimates restore m1
.                         mchange_mi std_same_prop_Sex "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi std_same_prop_AgeGrp4 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi std_same_prop_Race5 "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi std_same_prop_BornUSA "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi std_same_prop_MarStat5 "continuous" "mi"
.                 {c )-}
.         
.                 if (`model_version' == 2 | `model_version' == 4){c -(}
.                         estimates restore m1
.                         mchange_mi UnEmpl "categorical" "mi"
.                         estimates restore m1
.                         mchange_mi PhysProb "categorical" "mi"
.                 {c )-}
.         
.                 if (`model_version' == 2){c -(}
.                         estimates restore m1
.                         mchange_mi RAT_UnEmpl "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi RAT_PhysProb "continuous" "mi"
.                 {c )-}
.                 if (`model_version' == 4){c -(}
.                         estimates restore m1
.                         mchange_mi std_same_prop_UnEmpl "continuous" "mi"
.                         estimates restore m1
.                         mchange_mi std_same_prop_PhysProb "continuous" "mi"
.                 {c )-}               
.         {c )-}
. {c )-}
{txt}
{com}. 
. if ("`model'" == "logit"){c -(}
.         svy: logit Suic i.St `model_eq', or 
.         estimates store m1 
. 
.         if (`model_version' <= 4){c -(}
.                 mchange `margin_demographics', amount(all) delta(100) statistics(all) decimals(7)
.         {c )-}
. 
. {c )-}
{txt}
{com}. 
. if (`model_version' == 5 | `model_version' == 6) {c -(}
.         estimates restore m1
.         margin_interact std_same_prop_Sex Female 5 `model'
.         estimates restore m1
.         margin_interact std_same_prop_AgeGrp4 AgeGrp4 5 `model'
.         estimates restore m1
.         margin_interact std_same_prop_Race5 Race5 5 `model'
.         estimates restore m1
.         margin_interact std_same_prop_BornUSA BornUSA 5 `model'
.         estimates restore m1
.         margin_interact std_same_prop_MarStat5 MarStat5 5 `model'
. 
.         if (`model_version' == 6) {c -(}
.                 estimates restore m1
.                 margin_interact std_same_prop_UnEmpl UnEmpl 5 `model'
.                 estimates restore m1
.                 margin_interact std_same_prop_PhysProb PhysProb 5 `model'
.         {c )-}
. 
. {c )-}
{txt}
{com}. 
. 
. log close 
      {txt}name:  {res}<unnamed>
       {txt}log:  {res}/N/project/suicide_study/pnas_replication/results/log/mi_1.smcl
  {txt}log type:  {res}smcl
 {txt}closed on:  {res}21 Aug 2020, 03:00:58
{txt}{.-}
{smcl}
{txt}{sf}{ul off}